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Abstract 
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avoiding strictly to introduce any phenomenological interaction potential of nonperturbative origin) , 
may be improved by allowing for a merely numerical solution of the corresponding eigenvalue problem. 
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Recently, Pineda and Yndurain jlj] presented a re-analysis of heavy quarkonia. Their investigation 
is based on the main assumption that bound systems of heavy quarks may be reasonably described by 
nonrelativistic kinematics and only the perturbative contribution to the quark-antiquark interaction 
potential V if all nonperturbative effects are taken into account by some appropriate correction to the 
energy. In order to describe a system of a heavy quark and antiquark, both with constituent mass m, 
forming a bound state with total spin s = or s = 1, Pineda and Yndurain consider the Hamiltonian 
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The perturbative contribution to the static quark-antiquark interaction potential, V P (r), is known up 
to and including the two-loop level J2j : 
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Here, the following notations have been adopted: denotes the strong fine-structure constant in 
the modified minimal-subtraction (MS) renormalization scheme. For a non-Abelian gauge theory for 
ri{ Dirac fermions, invariant w. r. t. gauge transformations forming a Lie group SU(iV) describing N 
colour degrees of freedom, the quadratic Casimir invariants read, for the fundamental representation, 

N 2 -l 
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and, for the adjoint representation, 

C A = N , 

if the generators of the Lie group SU(iV) are normalized such that the second-order Dynkin index of 
the fundamental representation is 



The dependence of the effective (running) fine-structure constant a(/j.) on the renormalization scale y, 
is described in terms of the Gell-Mann-Low (3 function according to 
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involving the well-known expressions for the (gauge-invariant) one-, two-, and three-loop expansion 
coefficients in the MS scheme H 
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The resulting dependence of the fine-structure constant a(/j.) on the chosen renormalization scale y, 
expressed in terms of the (standard) scale parameter A, reads up to and including the three-loop level 
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7e is known as Euler-Mascheroni constant. In the case of quantum chromodynamics, clearly, N = 3. 
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Now, in order to stick to an entirely analytical analysis and following the philosophy developed in 
an earlier treatment [f| of heavy quarkonia, in Ref. Q the static potential Vq(t) ^ s split, according to 
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with the effective fine-structure constant 
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and an obvious remainder involving logarithms of the radial coordinate r, 
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The eigenvalue problem for the Coulombic Hamiltonian 

H=2m- — A + V(r) 
m 



is solved exactly, yielding, for instance, for the ground state, the energy eigenvalue Eq 
with the Coulomb binding energy 
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The "non-Coulombic" part H—H of the Hamiltonian (|l|) is treated perturbatively. Counting carefully 
the powers of a yields analytical expressions for the eigenvalues of H correct up to and including the 
order a 4 . Nonperturbative effects are incorporated by adding the Leutwyler-Voloshin correction S, 
which involves the gluon condensate (a G 2 ) — 0.06 ±0.02 GeV 4 . For the ground state, this correction 
amounts to the energy shift 

624 7T (a G 2 ) m 
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By inversion of the expression for the (ground-state) bound-state mass emerging from this procedure, 
the corresponding quark pole mass m is computed. For instance, with A(n,{ = 4) = 0.23 t° ' l GeV and 
/U = V&.632 ± 25 % GeV, implying a(/z) = 0.246 and a(/i) = 0.386, the experimental value || of the T 
mass, M(T) cxp = 9.46037± 0.00021 GeV, translates into the b quark mass m b = 5.001 ±g;^ GeV §. 

However, a perturbative treatment as implied by the splitting Q) of the static potential ^ P ( r ) ^ s 
by no means mandatory, obligatory, or even desirable. We may also adopt the following point of view. 
Given the operator H defined by Eq. (^) (accurate up to a certain order in a), compute (numerically, if 
necessary) its discrete spectrum, i. e., the set of eigenvalues, irrespective of the involved powers of a. 
Of course, the terms in Eq. (|l]) proportional to 

-A 2 , +-A, 
r 

and, if the effective coupling strength multiplying this term exceeds some critical value, also the term 
in Eq. (|l]) proportional to 

1 

render the operator H unbounded from below and have therefore to be treated perturbatively anyway. 
The Hamiltonian 



H = 2m- —A + V^(r) 
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on the other hand, may certainly be analyzed without adhering to some perturbative approximation. 
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In order to obtain a first idea of the differences brought about by these two approaches, let us start 
by considering only the Hamiltonian (|3|). The perturbative calculation is straightforward. Introducing, 
for notational brevity, the generalized Bohr radius 



a(/i) = 



the expectation values of the non-Coulombic interaction V(r) w. r. t. the ground state of H (indicated 
by the subscript 0) may be evaluated with the help of the relations (see also Appendix B of Ref. B) 
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The nonperturbative evaluation of H is performed with some numerical procedure^] developed for the 
treatment of the nonrelativistic Schrodinger equation [Q . In this way, we find, for the parameter values 
used in the second of Refs. |jj and focusing our interest to the ground state, for the Coulomb binding 
energy, e = —0.33146 GeV, for the expectation value of V(r), (V{r))Q = — 0.I37I4 GeV, and thus, 
for the perturbatively calculated ground-state energy,^ E p = 2 m + eo + {V(t))q — 9.5334 GeV, while 
the numerically computed "exact" lowest eigenvalue of the operator (||) is Eq P = {H)o = 9.5198 GeV. 
Consequently, for the parameters of Ref. |l| the difference in the lowest bound-state mass predicted by 
the Hamiltonian within perturbative and nonperturbative approaches is Eq — Eq P = 13.6 MeV. 
Unfortunately, this discrepancy is roughly 65 times larger than the experimental error on the T mass. 
Phrased the other way round, the perturbative ground-state eigenvalue E p of H can be reproduced by 
the nonperturbative evaluation of H for a mass m of the bound-state constituents of m = 5.008 GeV. 
For the b quark mass of Ref. |jj , the theoretical error different from the one induced by variation of the 
renormalization scale /i, attributed to neglected higher-order perturbative as well as nonperturbative 
corrections, is estimated to be ±0.006 GeV. Obviously, this error is entirely consumed already by the 
difference of the masses of the bound-state constituents obtained by perturbative and nonperturbative 
evaluations of H. 

In view of these findings, let's try to improve the theoretical value of the b quark pole mass nib by 
approaching the part H of the Hamiltonian (Q) nonperturbatively. The numerical computation of the 
expectation values of the operator H—H is considerably facilitated by the following two observations: 

• The eigenvalue equation for the "toy-model" (Hamiltonian) operator H defined in Eq. (Q) reads, 
for some generic (energy) eigenvalue E and its corresponding eigenstate |V>) of H, H\\jj) = E \ip). 
Thus the expectation values of all terms in the Hamiltonian (Q) involving the Laplacian A, taken 

w. r. t. \ip), may be evaluated by substituting A according to A|V>) = m 2 m + V p (r) — E 

• The expectation value of the S function entering in the "spin-spin term" of the Hamiltonian (|f|), 
taken w. r. t. \tp), is the modulus squared of the corresponding wave function ip(pc) at the origin: 
(ip\S^ (x)|V>) = IV^O)! 2 . For states with vanishing orbital angular momentum I (the so-called 
"S waves"), |V>(0)| 2 may be expressed in terms of the first derivative of the relevant interaction 
potential V(r) w. r. t. the radial coordinate r according to (for a derivation, see, e. g., Ref. 0) 
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In this way, the lowest eigenvalue of the Hamiltonian H may be computed from the expression (where 
the subscript of the expectation values indicates, as before, the ground state of the Hamiltonian H) 
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1 The desired accuracy of the (numerically determined) bound-state energies and wave functions may be adjusted in 
the routine used for the solution of the Schrodinger equation. For the present analysis, the uncertainty of these energies 
has been required to be less than 10 — 7 GeV = 100 eV. This accuracy should be, by far, sufficient for our purposes. 

2 For the perturbative treatment of H , we truncate the Rayleigh-Schrodinger series for Eq at lowest non-trivial order 
in V(r). Inclusion of the next order [|lj reduces the observed discrepancies but docs not change qualitatively our findings. 
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Adding the nonperturbative shift SEq gives our final result for the ground-state energy: £q = Eq+SEq. 
For the parameter values of Ref. |jj and a b quark mass of mb = 5.001 GeV, this expression entails the 
bound-state energy £o = 9.4953 GeV. Hence, the error of the predicted T mass brought about by the 
perturbative approximation in the treatment of the Hamiltonian ([l]), which amounts to 35 MeV, is of 
the order of the hyperfme splittings in the bottomonium system. For the latter, the second of Refs. |l| 
quotes, e. g., M(T) - M{-q h ) = 46.6 ±\$* 7 MeV. 

Consequently, a re-evaluation of the b quark mass appears to be in order. Fitting the ground-state 
bound-state energy £q, for the numerical values of the parameters A, fi, and (a G 2 ) adopted in Ref. jl[], 
to the experimental mass of the T, we obtain for the pole mass of the b quark 

m b = 4.983 GeV . 

The errors caused by the uncertainties of A and (a G 2 ) and by the variation of /i should be practically 
the same as in Ref. Q. Hence, dropping the requirement of analytical accessibility of the Hamiltonian 
(Q) reduces the extracted b quark mass by some 18 MeV. A very similar result is expected to be found 
for the determination of the c quark mass. We arrive at the conclusion that in Refs. pj the theoretical 
errors on the quark masses have been somewhat underestimated. 

Note that the above considerations apply, of course, also to the analysis presented in Ref. Q . Note 
also that the correct numerical factor in the numerator of the first term in the expression for the 
three-loop /3 function coefficient /?2 differs slightly from the one used in Refs. [Q. 
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